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ABSTRACT 

Ha emission from supernova remnants (SNRs) implies the existence of neutral 
hydrogens in the ambient medium. In the precursor of an SNR shock modified 
by cosmic rays (CRs), upstream plasmas are pushed by the CR pressure, but 
neutral particles are not, so that the relative velocity appears and some neutral 
particles become pickup ions by the charge exchange process in the precursor. 
We investigate how the pickup ions generated in the precursor affect the shock 
structure and the particle acceleration. If the CR pressure is larger than 20% of 
the shock ram pressure, the compression of the subshock becomes smaller than 
that without pickup ions because of the pressure of the pickup ions. Moreover, 
even if the shock is modified by CRs, the total compression ratio can be smaller 
than 4. In addition, the pickup ions play an important role for the injection into 
the particle acceleration. If the shock is a quasi-perpendicular shock and if the 
the multiply reflected ion acceleration occurs, the CR spectrum can be harder 
than that of the test particle diffusive shock acceleration below GeV. 

Subject headings: acceleration of particles — cosmic rays — plasmas — shock 
waves — supernova remnants 



1. Introduction 

Supernova remnants (SNRs) are observed in various bands from radio to gamma-ray. 
X-ray observations provi de the evidence that electrons are accelerated to highly relativistic 



energies in SNR shocks (iKoyama et al.lll995l ). Moreover, it is believed that ions are also 



accelerated to the "knee" (~ 10 15 eV) of the Galactic cosmic ray (CR) spectrum providing 
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most of the Gal actic CRs. The most popular acceleration mechanism is diffusive shock acce l- 
eration (DSA) flAxford et al.lll977l : iKrvmskv Ill977i : Iiellll978l : blandford fc Ostriker! \l97^ \ . 
To make the Galactic CRs, on average a few tens of percent of the explosion energy of a 
supernova has to be converted to the CR energy so that accelerated C Rs change the shock 
structure to make the precursor in the upstream (IDrury fc Volklll98ll ). On the other hand, 
the existence of neutral particles has been identified in many young SNRs from the detectio n 
of Ha emissions. A recent review of Ha emission from SNRs can be found in iHengl ( 120101 ). 
The neutral fraction of the interstellar matter around SNRs is often found to be order of 



unity. For example, iGhavamian et al.l ( 120001 . 120021 ) have shown that the neutral fraction is 
about 0.8 and 0.1 in Tycho and SN1006 environments, respectively. Since SNR shocks are 
collisionless, neutral particles can penetrate the shock front freely without deceleration. As 
a result, a significant relative velocity exists between the plasma and neutral particles in the 
shock downstream region when the shock is unmodified by accelerated CRs. Through the 
charge exchange process in the downstream re gion, hot protons are conve r ted to hot neu 



tral h ydrogens which emit broad Ha emission (j Chevalier fc Raymond! Il978l ; IChevalier et al. 



19801 ). Conversely, the cold neutrals are converted to cold protons which are not equilibrated 
with the downstream thermal protons due to the Coulomb collision because the relaxation 
time scale t T 



t PP ~ 



3 x 10 5 yr 



Vrel 



.lcm- 3 / V 3000 km s" 1 , 
is much larger than the age of young SNR (~ 10 3 yr), where n p and t> re i are the proton 
number densi ty and the relative y elocity between the cold protons and the thermal protons, 
respectively. lOhira et al.l (j2009bl ) have shown that the charge-exchanged ions amplify the 
magnetic field and change the shock structure in the downstream. 

In this letter, we propose that charge-exchanged ions are also important in the precursor 
of the shock modified by the CR because the relative velocity between hydrogens and protons 
exists in the precursor. Although the neutral hydrogen damps the plasma wave indispens- 
able for DSA in the upstream region, the wave amp litude is still kept to be strong enoug h 
through the excitation by CRs near the shock front (jDrury et al. 1996; iReyille et al. 2007 ). 



Henc e, CRs can be accele r ated i n the partially ionized gas ( IHelder et al 



2009 



Lee et al. 



2010l ). In fact, iLee et al.l (120101 ) have observed Ha from the precursor in the SNR Ty 



cho. The charge-ex changed ions excite plasma waves (IWu fc Davidson 



(Wu & Davidson 


1972; 


Lee & Ip 


1987; 


(Williams & Zank 


199^ 


: ). As the re- 



suit, the plasma waves amplify the upstream magnetic field and the charge-exchanged ions 
become isotropic pickup ions. To estimate the efficiency of the charge exchange in the pre- 
cursor, we compare the precursor length £ pre with the mean free path of the charge exchange 
^c.e. = (%>0"c.e.)~\ where oc.e. is the cross section of charge exchange which is almost con- 
stant o"c.e. ~ 10~ 15 cm 2 for t> re i < 3000 km s _1 and a steeply decreasing function of t> re i for 
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v TC \ > 3000 km s -1 ( Heng fc McCrayl 120071 ). We here estimate the precursor length as the 
diffusion length of CRs £d = D(p)/uo, where p, D{p) and u are the CR momentum, the 
diffusion coefficient of CR and the shock velocity, respectively. Assuming the Bohm diffusion 
for CRs D{p) = c 2 p/3eB and cxc.e. = 10~ 15 cm 2 , we obtain 

T^T = 3 X 10_2 (^) (t^) (3OOO km s- 1 ) (37Ig) ' (2) 

where B is the magnetic field strength. Thus, we expect that more than a few percent of 
hydrogens are converted to the pickup ions in the precursor. Because on average 10 -4-5 
fracti on of the upstrea m protons are needed to be injected into DSA to supply the Galactic 



CRs (IVolk et al.ll2002l ). the pickup ions are clearly important for the injection into CRs and 



for the modification of the shock structure. 

We first provide a simple model of the modified shock structure (Section 2). We then 
solve the distribution function of the pickup ions (Section 3) and discuss impacts on the 
jump condition of subshock (Section 4) and the particle acceleration (Section 5). 



2. A simplified model of the modified shock structure in the precursor region 



To calculate the distribution of pickup ions, we need to know the relative velocity 
between protons and hydrogens in the precursor. In this section, we make a simplified 
model of the precursor structure of modified shock. We assume that the pickup ions do not 
change the precursor structure, that is, the number of pickup ions is smaller than that of 
thermal protons. The shock structure modified by the CR pressure has been calculated by 



various methods flDrurv fc Volklll98ll : Ellison et al.lll99CH : iBerezhko fc Volklll997l : lBlasill2002 



Kang et al.ll2002l ). In order to retain an essential feature of the structure of modified shocks, 
we determine the approximate structure of the precursor. Mass and momentum conservation 
laws at the shock rest frame are 



Pqu = p(x)u(x) 
PquI + P gi0 = p(x)u(x) 2 + P g (x) + P C r(x) 



(3) 



where p, u, P g and Pgr are the plasma density, flow velocity, the plasma pressure and the 
CR pressure, respectively. Here, subscript means quantities far upstream and no subscript 
means those in the precursor, and x is the coordinate of the direction along the shock normal 
(x = and x = —00 are the position of the subshock and far upstream, respectively). 
Neglecting the gas pressure, the velocity in the precursor is 



Uo 



Poul 



(4) 
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where Pgr is given as 

/"Pmax 1 

Pcr(x) = / ^PcfcRiP, x)4iTp 2 dp , (5) 

J mc " 

where /cr(p, is the distribution function of CR, and we assume that the pressure of 
non-relativistic CR is negligible. To make the expres sion simple, we use fcn(p,%) as the 



approximate solution of test particle DSA (IBlasil 120021 ) 



/cr(p, x) oc p s exp ( ^rr j ~ P s ©(4(p) + a?) , (6) 

where 9 is the Heaviside step function. From Equations 01]), dSJ) and flBJ), we obtain the 
velocity structure of the plasma in the precursor 

71^ f Tmag-(k'|/^d,mc) 4 ~ S In I I ^ \ 

W = j _ t x J l«d,max ^ |X| > «d,mcj , ? v 

U ^ \l (N<4,mc) 

where 4,mc = 4(wc), £ = Pcr(0)/po^O' an d 7max is the maximum Lorentz factor of CRs. 
We presume that £ is order of 0.1 to be consistent with the general scheme of CR modified 
shock and the Galactic CR production. Therefore, the velocity deceleration in the precursor 
region is around u — u(0) = £u ~ 300 km s -1 . For s < 4, the precursor scale £ pm is 
4,max = 4(Pmax) and for s > 4, 



3. Distribution function of pickup ions 

In this section, we calculate the distribution function of the pickup ions at the subshock, 
f(v,x = 0). Charge-exchanged ions excite plasma waves and make a bispherical distribution 
by the pitch angle scattering with the excited waves. Here, we approximate the bispherical 
distribution as an isotropic distribution at plasma rest frame since the Alfven Mach number 
(M A ) of SNR shocks is high. 

From equation (J2J), £ pre is usually smaller than £q.e.- Then, pickup ions are produced in 
the whole precursor region and the flow velocity of hydrogens is uq in the precursor. Hence, 
the relative flow velocity between the hydrogens and protons Au is 

Au(x) = u — u(x) . (8) 

We here assume that the thermal velocity of protons and hydrogens is smaller than Am. 
The validity of this assumption is discussed later. The transit time of pickup ions over the 
precursor region is r = £ prc / (1 — £)w which is less than £c.e./ (1 — u 0i so that 

T < 01 yr (d^r (y ) ' ( 3000I,- )"' • (9 > 
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From equation (JTJ, the equilibration time of the Coulomb collision in the precursor is 



The recombination time t rec is 

= lO^yr^j ^ | , (li; 

where we use the electron temperature T e because the electron thermal velocity is comparable 
to or larger than Am. t pp and t rec are much larger than the transit time, so that we neglect 
the collision and the recombination processes. Moreover, we neglect the charge exchange 
process between pickup ions and hydrogens because of £ pre < £q.e.- Hence, we consider only 
the charge exchange process between hydrogens and thermal protons. 

Then, the steady-state Boltzmann equation for pickup ions in the shock rest frame is 

, s df / \ / \ a / \S(v — Au(x)) , . 

u{x)— = n R (x)n p {x)ac.EAu{x) , (12) 

where v is the velocity of pickup ions at plasma rest frame and tzh is the hydrogen number 
density. The solution is 

f(v, x = 0) = -r^-^- f dxn H (x)^-S(v - Au(x)) , (13) 

where £q.e. = (%>,o0"c.e.) _1 and n Pi o is the proton number density far upstream. Then the 
distribution function of the pickup ions at the shock front is 

f{v,x = 0) = f 1 (v) + f 2 (v) , (14) 

where fi and f 2 are the distribution functions of the pickup ions ionized in the x < —£d,mc 
and x > —idmc, respectively, and they are 



»H 4,mc £ 7max ~ 1 ( J_\ [ is _ ( 4-s ^ ( V 

^v 2 £ c .e. (1 - v/u ) 2 £u (4 - s) \^u J \ 7max l7max ) \£u 



3-s 
" 4-s 



£d,mc £ 

4irv 2 £c.e. (i - 

where we approximate as a constant. For 7 max ^> 1, fi is approximately expressed by 



f^) = ^7^JT ± TvA--^o) , (16) 



f() ~ Wh ^ d ' max ^ 1 ( _v_ \ 

h[V) ~ Anv 2 £ C . E . (l-v/u ) 2 ^u \A-s\ \£u ) * 



— 3— s 

7^{l + 7^(|^)}~ 4 " S («>4) 



3- 

" 4-s 



(17) 
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4. Effect of pickup ions on the jump condition of the subshock 



In this section we discuss the jump condition of the subshock in the presence of pickup 
ions. As is well known for the termination shock of the solar wind (up ~ 400km s" 1 ) , large 



number of pickup ions change the jump condition of shocks ([Richardson et al.l 120081 ). The 
high pressure of pickup ions makes the Mach number and the compression ratio small. In 
addition, there is no guarantee that the behavior of the pickup ions in the sh ock is the same 
as that of gas with adiabatic index 5/3 (IFahr fc Chalovl 120081 ; IWu et al.ll2009l ). Interestingly, 
the peculiar behavior (adiabatic index > 5/3) of pickup ions may allow the possibility that 
the total compression ratio can be smaller than 4 and the spectrum of CR can be softer than 
predicted by the test particle DSA, even if the shock is modified by the CRs. The pressure 
of pickup ions at the subshock front is 



pi 



dw 2 f(v) = Ppi,! + Ppi, 2 



18) 



where Ppi,i and Ppi,2 are partial pressures due to fi and fi-, respectively. Approximating the 
third term of Equation f[TT|) £/(! — v/u ) 2 by £/(! — £) 2 , we obtain 



pi,i 



pu z 



1 Ph 
3 po 



ic.E. 



(1 



_3 I 3 

5-s ^ 9-2s 



,mc -,13 — 3s 

in 



1 



13-3s 
3 



£ CE /max 

^d,mc 1 



13-3s 



— h 

9-2s ~ 5-s 



) (s < 4) 

- 1) (4 < s < 13/3) 
(s > 13/3) 



(19) 



<C.E. 3a- 13 

where pn is the density of hydrogens. From Equations ffl6l) and (fT8]) . we obtain 

Ppi,2 1 PH ^d,mc / C X 



3 po 



-C.E. 



1"£ 



(20) 



If the modification of the shock is not strong enough (^ < 0.2), the pressure of pickup ions 
is negligible (Ppi/pu 2 < 0.01) , so that the effect of pickup ions on the precursor structure 
is negligible and our treatment in section [2] is valid. If s < 13/3, pu ~ p , £ ~ 0.5 and if 
most hydrogens become the pickup ions in the precursor, from equation f fl9|) and ( 120]) . the 
pressure of pickup ions account for about a few tens of percent of the ram pressure at the 
subshock, that is, the Mach number at subshock is almost 1. Therefore, the compression 
ratio at the subshock is almost 1, so that the total compression ratio can be smaller than 4. 
Although our treatment in section |2] is not valid in this case, we qualitatively can expect the 
reduction of subshock compression and the strong magnetic field amplification. 

To calculate the jump condition of the subshock rigoro usly, we have to study t he behav- 



ior o f the pickup ions when they pass over the subshock (IMatsukiyo et al.ll2007t IWu et al. 



20091 ). Although it should be investigated by particle in cell simulations, possible effects are 
mentioned in the next section. 
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5. Effect of pickup ions on particle accelerations 



As already discussed, the total compression ratio may become smaller than 4. Therefore, 
even when the SNR is young, that is, the Mach number of the shock is large, the spectral 
index of CRs may be softer than 2. 

At present, the injection process for CR accelerations at shocks is still an open problem. 
It is regarded that some particles above the threshold velocity are injected into the accel- 
eration process. As discussed above, the heating of plasma at the subshock is weak when 
the pickup ions drain a large fraction of the shock kinetic energy. Therefore, the pickup 
ions have the velocity much larger than that of thermal plasma and are important for the 
injection into CRs. This is a similar process to the solar wind termination shock where the 
pickup ions have been considered as the origin of the anomalous CR. 

The behavior of the pickup ions at the subshock of SNR has rarely been investigated. 
While the adiabatic accelerati on is an important process, in thi s letter, we discuss the multi- 
ply reflected ion acceleration (IZank et al.lll996t iLee et al.lll996l ) as another possible process. 



When the shock is the collisionless quasi-perp endicu 



tential to reflect some upstream thermal ions (ILerov 



ar sh ock, there is an electrostatic po- 



1983). The reflected thermal ions can 



penetrate the shock front after the reflection. iLeroyl (119831 ) estimated the amplitude of 
the electrostatic potential in the steady-s tate framework. However, for the high Ma shock, 
the shock is not stationary (IQuestl Il985l ) and the amplitude of time-averaged electrostatic 
potential is estimated by 

e<p = r]-mu 2 , (21) 



where rj is an order of unity parameter (IShimada fc Hoshindl2005l ). In the rest frame of the 
precursor plasma, particles satisfying v x < — A£uo are reflected by the electrostatic potential, 
where v x is the velocity of the shock normal direction and A = £ -1 (l — £)(1 — r/ 1 ^ 2 ). There 
are many pickup ions satisfying the reflection condition compared with thermal ions, and 
the pickup ions can be reflected by the electrostatic potential many times. The reflected ions 
drift to the direction of the motional electric field Uo x B. As a result, the reflected ions are 
accel erated. This acceleration mechanism is so-called shock surfing acceleration ( jSagdeev 
19661 ). The number of the pickup ions reflected by the potential n re f is 



'•Tef.s 



/■ftio p6 

2tt / dvv 2 fi{v) / dd'smd' 

J A£u JO 



(22) 



where 9 is defined by cos# = A^uq/v, and % = 1 and 2. Since the total number density of 
the pickup ions npi ti = 2n rcii (A = 0), from Equations (fT6l) . (jPTl) . and (1221) . the reflection 
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coefficient of the pickup ions Ri = n re f,i/?"ipi,i are 



Ri 
R2 



X 



l-A(5-s)-(s-4,)(l+A^)' 



/max 



l-(*-4)7S 



(a<4) 
(a>4) 



(23) 
(24) 



Figured] shows the reflection coefficients where the red, green and blue lines show Ri(s = 3.5), 
Ri(s = 4.5, 7 max = 10 5 ) and R2, respectively, the solid, dashed and dotted lines show 77 = 0.9, 
7] = 0.5 and 77 = 0.1, respectively. There exists the critical £ at which the reflection coefficient 
becomes 0. From the condition A — 1, this critical £ cr is 



1 



V 



1/2 



V 



1/2 



(25) 



When £ > £ cr , there are reflected pickup ions, so that it is possible that the multiply re- 
flected ion acceleration occurs. When the shock normal component of the Lorentz force 
becomes larger than the electric force of the electrostatic potential, the particle escapes to 
the downstream. So that the maximum velocity of accelerated particles is given by 



e B 1 ~ e- 



v ramp 



(26) 



where L ramp is the thickness of the electrostatic potential and B± is the magnetic field 
perpendicular to the shock normal direct ion. The observatio n of the termination shock 
of the solar wind shows L ramp ~ c/a; PiP (IBurlaga et al.l 120081 ). From Equation (126]) . the 
maximum velocity is 



77 /uy u PjP 



2 \c) O 



c.p 



2[? 



u 



3000 km s" 



1 cm 



1/2 / B 



3 fiG 



-1 



(27) 



Therefore, the multiply reflected ion acceleration at SNR shocks can make relativistic CRs 
and is important for the injection mechanism. This acce leration makes a harder spectrum 
than that of DSA (IZank et al.lll996l ; iLipatov fc Zanklll999l ). so that this m ay explain that th e 
observed 7-ray spectrum from W51C is harder than expected for DSA (jAbdo et al.l 120091 ). 
However, since the high Ma (Ma > 100) collisionless shocks are not well understood, i.e., 
the shock structure may be different from the solar wind termination shock, it should be 
investigated whether the multiply reflected ion acceleration really occurs or not at SNR 
shocks. 
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Discussion 



We discuss some important processes neglected in the above investigations. If consider- 
able hydrogens are converted to the pickup ions as discussed in section [31 they may change 
the structure of the precursor because of their large momentum flux. This issue will be 
addressed in future work. 

In Section |3l we neglected any processes other than the char ge exchange. If the upstream 



plasma is strongly heated by unstable waves in the precursor ( 


Bell 


2004; 


Ghavamian et al. 


2007; 


Rakowski et al. 


2008; 


Niemiec et al. 


2008; 


Wagner et al. 


2009 


■ 


Ohira et al. 


2009a 


), that 



is, if the temperature becomes T > m(£w ) 2 , neutral hydrogens are ionized by the electron 
scattering and the ionized particles can not be distinguished from thermal ions because 
thermal ions have the velocity comparable to that of the pickup ions. iRakowski et al.l (120091 ) 
showed that the collisional ionization occurs in the upstream of SNR. 

In Section El we considered only the case £ pTe < Iq.e.- If the magnetic field is not so 
strong and s < 4, the precursor scale is the diffusion length of the maximum energy CR 
and £ pre > Ic.r (see equation (El) and (J7D). For £ pre > £c.e., the multiple charge exchange is 



important (Ivan Adelsberg et al.ll2008l ). To take account of it and collisional ionization, we 
must fully take into account the velocity dependence of the various atomic cross sections. 
We will make full calculation in future work to investigate these issues. 



Raymond et al.l (120081 ) predicted that the pickup ions introduce a non-Gaussian broad 
Ha line from downstream region. We also predict it from upstream region. So far, in 
calculations of Ha from SNR, one has treated the pickup ions and thermal plasma as one fluid, 
and assumed that the pickup ions do not affect the plasma flow. As mentioned previously, 
the pickup ions should be treated as a separate component. 



7. Summary 

In this letter, we have investigated the effects of neutral particles on the shock modified 
by the CR at SNRs. More than a few percent of hydrogens are converted to pickup ions by 
the charge exchange in the precursor. In Section [31 assuming that the pickup ions become 
isotropic in the rest frame of the precursor plasma, we calculated the expected distribution 
of pickup ions at the precursor. The pickup ions amplify the magnetic field and heat the 
upstream plasma. Moreover, the pickup ions have a large pressure in the precursor and 
change the jump condition of the subshock when the upstream plasma is strongly pushed 
by the CR pressure (£ = -Pcr(0)/poWq > 0.2). In addition, because of the large velocity 
of the pickup ions, they are preferentially accelerated by DSA. If the shock is the quasi- 
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perpendicular shock and if the multiply reflected ion acceleration occurs, it is possible to 
make harder spectrum of CR than predicted by the test particle DSA below GeV. Hence, 
the pickup ions are important for the injection into CR. 

We thank T. Terasawa and M. Scholer for useful comments about pickup ions. We also 
thank the anonymous referee for valuable comments to improve the paper. This work is 
partly supported by Scientific Research Grants (O.Y.: 21684014, and F.T.: 18542390 and 
20540231) by the Ministry of Education, Culture, Sports, Science and Technology of Japan. 
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Fig. 1. — The reflection coefficient of the pickup ions. The red, green and blue lines show 
Ri(s = 3.5), R\(s = 4.5,7 max = 10 5 ) and R 2 , respectively. The solid, dashed and dotted 
lines show rj = 0.9, rj = 0.5 and rj = 0.1, respectively. 



